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Equations for the radiation field of an arbitrary distribution of source currents embedded 
in a lossless magneto-ionic medium and radiating at a frequency that is high relative to the 
electron gyrofrequency of the medium are derived using the QL approximation. The field 
is given directly in terms of the components of radiation vectors determined by an integra- 
tion over the source currents. The application of the equations to arrays of similar sources 
is discussed. 



1. Introduction 

The subject of radiation from current sources in anisotropic media has received a great 
deal of attention in recent years. One reason for some of this activity is the current interest 
in space technology, plasma physics, and the problem of transmitting through the reentry 
plasma sheath [Bachynski, 1961]. 

Bunkin [1957] has published a formal solution for the radiation field of monochromatic 
current sources in an infinite homogeneous linear gyroelectric medium. In an earlier report 
Abraham [1953] attacked the problem of radiation from sources in a magneto-ionic medium 
by a somewhat different method. Since the publication of Bunkin's paper, a large number of 
papers related to this subject have appeared, including those of Chow [1962], Ford [1961], 
Kogelnik [1959], Kuehl [1962], and Shore and Meltz [1962]. 

The formal solution to the problem of electromagnetic radiation in gyroelectric media as 
given by Bunkin [1957] is incomplete in so far as it is necessary in order to obtain the field 
vectors, first to solve a transcendental equation [Bunkin, eq 3.13] and, second, to calculate the 
components of a matrix using this solution. While it does not appear to be feasible to obtain 
a general solution to the transcendental equation in question, solutions are possible in certain 
special cases. Bunkin, in the paper referenced, gave some results for the field of a dipole in 
two special directions, parallel to the static magnetic field of the medium and perpendicular to 
it. Kuehl [1962] lias obtained solutions in several cases where special forms of the dielectric 
tensor are assumed. In one of these cases the limiting form at high frequency of the dielectric 
tensor of a lossless magneto-ionic medium (cold magnetoplasma) was used. 

The subject of this paper is a case somewhat more general than the one for high frequency 
considered by Kuehl. While the same form of the dielectric tensor is assumed, and the rela- 
tive gyrofrequency is likewise taken to be small, the relative plasma frequency is restricted 
only by the stipulation that the source frequency exceeds the plasma frequency. Also, terms 
of the first order in the relative gyrofrequency are retained in the expression obtained for the 
electromagnetic field. This expression is given in a form analogous to that commonly used 
[Schelkunoff, 1939] for the calculation of radiation from current sources in isotropic media. 
Both electric and magnetic type sources are considered, and the application of the expression 
to arrays of similar sources is discussed. 



i Work done in part ial fulfillment of the requirements for the Ph. D. degree in Electrical Engineering at the University of Maryland. Dr. 
H. Reed acted as thesis advisor. 

2 Present address: IIT Research Institute, Electromagnetic Compatability Analysis Center, USNMEL, Annapolis, Md. 
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2. Formal Asymptotic Solution 

In this section the formal asymptotic solution for the field due to sources embedded in a 
lossless magneto-ionic medium is given. 

2.1. Maxwell's Equations 

In a magneto-ionic medium, Maxwell's equations for harmonic fields varying as exp 
(jut) are 

-jcoeeE+VXH=J (la) 



-VXE—ju f xH=M. 



(lb) 



A fictitious magnetic current density M is included [Schelkunoff, 1943] and rationalized mks 
units are used. E, H, and J" have their usual significance as the electric intensity, the magnetic 
intensity, and the electric current density respectively. The constants e and /x are the electric 
permittivity of free space and the magnetic permittivity of free space respectively. The 
relative electric permittivity, e, is a tensor. If Cartesian coordinates x, y, and z are used with 
the 2-axis directed along the static magnetic field in the medium, the tensor e takes the form 
[Suhl and Walker, 1954]: 

a —jb 0" 



Jb 







(2) 



In the case under consideration, the quantities a, b, and c are assumed to be given by the 
usual expressions in the magneto-ionic theory [Ratcliffe, 1959] with the collision frequency 
set equal to zero: 

a=l-X/(l-P) (3a) 



b=XY/(l-Y 2 ) 
c=\-X. 



(3b) 
(3c) 



where X is the square of the relative electron plasma frequency, X=cof/co 2 , and where Y is 
the relative gyrofrequency, Y=o) g /u. 




Figure 1. Components of radius vector R. 




Figure 2. The components of the vector N. 
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2.2. Method of Solution 



The components of e are independent of the coordinates in view of the assumed homo- 
geneity of the medium, although they do depend on w. It is possible to obtain an asymptotic 
solution to (1) by using the Fourier transform pair 



+ 00 

f(N)= f| f f(R)e im - R dxdydz 



(4a) 



+ O0 

f(R) =^3 fjfj(N)e-*».*d P dq ds 



(4b) 



in which R is the radius vector with Cartesian components x, y, and z as shown in figure 1, 
in which N is a vector whose Cartesian components are the real parameters p, q, and s as shown 
in figure 2, and in which N- R is the scalar product of these two vectors. The transform pair 
(4) is valid for any positive real value of the constant k. Here k is taken equal to coVm^ 

The solution of (1) is accomplished by applying the transform (4a) to both sides of (la) 
and of (lb) 3 to obtain a set of six algebraic equations in six unknowns 
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(5) 



where ^^VmA- By using partitioned matrices, eq (5) can be abbreviated as 



A 

-N 



A " 

N 




E 




—Jljut 


A 
/_ 




_vH_ 




-M/jk_ 



(6) 



where N is defined as 



N'= 



-s 
+ s 

L-2 +v 



+q.~ 
-v 





(7) 



and / is the identity matrix. 

The next step in the solution is to premultiply (6) by the inverse of the dimensionless matrix 



(8) 



This inverse can be written as 









- A 
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A A 






as 
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A 

-C(p,q, 





A 

D(p,q,s)_ 



(9) 



s The usual procedure of first eliminating the vector H from the pair of equations does not simplify the algebraic computations. 
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where the matrix on the right-hand side of (9) is the adjoint of the matrix (8) and A is the 

A 

determinant of the matrix (8). 4 The form of A has been given elsewhere [Bunkin, 1957, eq 

A A A 

5.3; Kogelnik, 1959, eq 21]. B, C, and D can be obtained by matrix multiplication [Marini, 

A AAA AA AAA A 

1963]: B=-AN, C=+NA, and D=NB+Al. 

The transform (4b) is then applied, giving formal integral expressions for E and H. The 
latter integrals can be evaluated by the method of residues and the method of steepest descent 
to give the radiation field. The procedure for this evaluation has been discussed elsewhere 
[Bunkin, 1957; Lighthill, 1960; Chow, 1962; Shore and Meltz, 1962] and will not be repeated 
here. The final expressions for the radiation field are given in the next section. 

2.3. Formal Solution for the Radiation Field 
In Cartesian coordinates the complete formal solution for the radiation field is 



E(R)=[jawA V +jkB W }g 



H(B)=[jkC Vo+j<*eD W ]go 



9 -jkN (R).R 

P 
-jkN (R)-R 



-[jccnA e V e +jkB e W e ]g e 



-[jW e V e +jueD e W e ]g e 



e -jkN e (R).R 
P 


(10a) 


e -jkN e (R).R 


(10b) 



where 



90,e~ J- 



±cos 6 

4x cos a 0> t 



and where 



n , e sin a . t 



sin 6 cos 2 (a 0>e — 0)[«' e — n 0ie ) cos (a 0fe — d)—2n' , 6 sin (a , e —0)} 

1 



[(c— a) 2 < e sin 4 a , e — 46 2 c« e sin 2 a , e — c)] 



Ao, e = A[p 0te , ^ 0i g, s 0f e \ 

A A 

Bo,e = B[po >e j ^o, ey s 0,e\ 
A A 

A A 

^0,e == -^[Po,ej Q_0,e) S 0,e\' 



(ID 

(12a) 
(12b) 
(12c) 
(12d) 



In eq (11) the lengths n (ao) and n e (a e ) of the vectors N (R) and N e (R) are obtained from the 
formula for the index of refraction 5 



<e(«) = 



2ac+(a 2 — b 2 — ac) sin 2 a±^(a 2 — b 2 — ac) 2 sin 4 a+46V cos 2 a 



2 (a sin 2 a+c cos 2 a) 



(13) 



In eq (12) the Cartesian components p , q 0j s and p e , q e , s e of the vectors N (R) and N e (R) are 
obtained from their spherical components 



Po, e=no,e(ao, e ) SU1 a 0> e COS ft), e 

<Zo, e=n>o, e («o, e) sin a Q , e sin Ot e 
r ,e=no,e(ao,e) sina , e 

So,e=%e(ao,e) COS O: , e . 



(14a) 
(14b) 
(14c) 
(14d) 



4 A is also equal to the determinant of the matrix e+A r 2 which premultiplies E when the usual procedure of first eliminating the vector //from 
eq (1) before applying the transform (4a) is followed. 

5 In eqs (11) and (13) the subscript corresponds to the upper sign while the e corresponds to the lower one. The abbreviations and e (for 
"ordinary" and "extraordinary") can be used here since the only case under consideration is the one where the quantities a, b, and c are given by 
eq (3) and a frequency higher than the plasma frequency is used [Ratcliffe, 1959, sec. 7.3]. 
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Vq, e and W 0> e are radiation vectors 



+™ 



Vv, e =fNe' m o,e™-*'J(R>)dx'dy'dz' (15) 

W , e = Nfe m o,eW'X'M(R')dx'dy'dz' (16) 

calculated by a direct integration over the current sources. 

The solution given in eqs (10-16) above depends entirely on the values assigned to (he 
azimuthal angles /3 and fi e , and to the polar angles a and a e . These are determined from 
the saddle point evaluation, which gives 

Po.e=<P (17) 

and 

7fc'(a ,e)Ma 0fg )=tan W e —6). (18) 

Equation (18) is the transcendental equation mentioned in the introduction upon which 
the entire solution depends. Kuehl [1962, eq (13)] has given explicitly the form taken by this 
equation when n^ e {a) given in (13) is substituted into it. If this form is specialized to the 
case of a lossless magneto-ionic medium by using the values given in eq (3) for the quantities'' 
a, b, and c, the eq (18) takes the form 7 

f AT1 ( _0W =FXF 2 sina 0iC COSQ:o, e 

iai U a 0)C v) 2 m 2 v(a , e )[l-Y 2 sm 2 a 0<c /2m 2 ±v(a 0<e )] ^ } 

In (19) above the abbreviations 

m 2 =l-X (20) 

K«)=VF 4 sin 4 a/4m 4 +F 2 cos 2 a (21) 

have been used. The quantity m may be interpreted as the length of the vectors N 0i e if the 
static magnetic field in the magneto-ionic medium is set equal to zero. 

If the values given in eq (3) are used in (13), then the Appleton-Hartree formula for the 
index of refraction in a lossless magneto-ionic medium results: 

nl e (a) = l-X/[l-Y 2 sin 2 a/2m 2 ±v(a)\. (22) 

3. Approximate Solution 

In this section the method of approximation used to obtain a solution to eq (19) will be 
indicated, and the resulting equations for the electromagnetic field will be given. 

3.1. Method of Approximation 

As mentioned previously, only the case where the wave frequency is high relative to the 
gyromagnetic frequency of the medium will be treated. Specifically it will be assumed that 
Fis small enough to permit the simplification of (19) to be carried out by expanding all quan- 
tities in powers of Y and dropping high order terms. 

In the case of an isotropic medium with Y equal to zero, the right-hand side of (19) vanishes 
and a 0tC can assume only the single value 0, which is the correct result. In the case under 



e The algebra is simplified if use is first made of the relation a 2 -b 2 —ac=a—c that exists between the quantities a, b, and c defined in (3). 
7 A factor o«/co appears to be missing from the corresponding equation given by Kuehl [1962, eq (46)]. 
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consideration, consequently, it is assumed that the angles a , e deviate only slightly from the 
limiting value d and that the deviation a , e — d, has an order of magnitude no larger than that 
of Y. In view of the expansion 

tan (a o ,e-0) = ao,e-e+(a o ,e-df/3+ . . . (23) 

eq (19) can be approximated by 

* -p XY 2 sin ao. e cos q , g ( ) 

a °' e ■*" 2m 2 <a ,e)[l-r 2 sin 2 a ,,/2m 2 ±^Ke)]' l j 

It remains to approximate sin a 0t e , cos a 0t e , and v(a 0<e ). The first two are approximated 
by writing a 0t e as a 0i e =0+ («o, e~ 0) and using appropriate trigonometric expansions. 

Some care must be exercised, however, when the quantity v(a) defined in (21) is approxi- 
mated. For most values of at, the approximation 

v(a)^Y\cos a\ (25) 

can be made, since the first term under the radical contains Y to the fourth power. If a is 
very close to ir/2, nevertheless, then the approximate value of (21) becomes 

v(a)^ Y sin 2 a/2m 2 +m 2 cot 2 a. (26) 

Equations (25) and (26) correspond to the QL and QT approximations respectively of the 
magneto-ionic theory [Ratcliffe, 1959, ch. 8]. In view of the assumption that Fis small, the QL 
approximation will hold over most of the region surrounding the sources. The QT approxi- 
mation is valid only over a small angular region in which rays from the sources are approxi- 
mately perpendicular to the static magnetic field. For this reason only the treatment of the 
QL case is included in this paper, although the QT case also can be handled in the same way 
[Marini, 1963]. 

3.2. Components of A^o, e (i?) 

The resulting approximate solution to (19) for the polar angle of N 0te (R) in terms of the 
polar angle of R is 

a ,e=OTXYsm B cos 0/2m 2 |cos 6\. (27) 

An additional term can easily be obtained in the expression (27) for a 0t e but the accuracy given 
is sufficient for present purposes. 

Using (27) above in (22), the length of N , e (R) is found to be given b}^ 



n%, 



= m 2 [l±^|cos^|-^(l+cos 2 ^)] (28) 



from which 



•o. e =m{l±g|cos0|-gJ[m*+(l-f)cos*0]|. (29) 



The azimuthal angle of N , e (R) is, as previously mentioned, 

&,.=*. (30) 

Equations (27), (29), and (30) above give approximate values for the spherical polar 
components of N (R) and N e (R) in terms of the spherical polar components of R. The Cartesian 
components are found from (14). 
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3.3. The Radiation Field 

The approximate equations for the radiation field can now be obtained by substituting 
the approximate values of the components of N 0te (R) found in the preceding section into eqs 
(10), (11), (12), (15), and (16). The components of the field vectors will be in Cartesian 
coordinates. 

AAA A 

It is evident that the matrices A 0t e , B , e , Co, o and D Q e are the Cartesian components of the 
tensors that relate the components of the radiation vectors to those of the field vectors. The 
use of Cartesian coordinates was dictated principally by the fact that the components of e can 
be made to assume the simple form (2) in that system of coordinates. Kuehl [1961] lias shown, 
however, that the equations for the electromagnetic field are more compact if coordinate sys- 
tems based on the planes of constant phase of the radiated waves are used. In the case under 
consideration, with a , e — d small, almost equal simplification is obtained by using spherical 
polar coordinates, and this will be done in what follows. 

One more simplification can be made. The absolute value sign that appears in (27) and 
(28) can be removed. This necessitates a change in notation, however, since the identity of 
the ordinary and extraordinary wave is switched in going from the upper half to the lower 
half space. 

With these changes, with the vector electric intensity written in spherical polar coor- 
dinates as 

"E D 



E(B) = 



E { 
IE J 



(31) 



and with a like notation used for the magnetic intensity, the radiation vectors, and the current 
densities, the equations for the field vectors are 



E(B) = (i/87rp)[co M Ai^ 
H(iJ0 = (i/8xp)[JfcC 1 V 1 +^^ 



(32a) 
(32b) 



where A 1>2 , B lt2 , C 1>2 , and D 1>2 are given by 

0; ±h sin 6; 

±hsin 8; -l±h(cos d— sin 2 B/X cos 8); 
_ jh sin 8; Tj+jh cos 6; 



Aj 2 — 



—jh sin 8 
dzj—jh cos 8 
■ 1 ± A (cos 0+sin 2 6/X cos 0)_ 



(33) 







0; jh sin 6; 


±h sin 8 








A 

Bi t2 =m 


— jh sin d; Tj; 
_ ±h sin 6; 1 T h sin 2 0/X cos 6 ; 

0; jh sin 6; 


-ITh sin 2 6/X cos 8 
Tj 

±h sin 8 








A 

C 12 =m 


—jhsind; ±j; 


l=FAsin 2 0/A r cos 8 






and 


" 0; 


_ ± h sin 6; -ITh sin 2 6/X cos 8 ; 
=F h sin 6 ; 


±3 

jh sin 8 




-. 


A 

Di, 2 = m 2 


TAsin 6; -l=F//(cos 6— sin 2 6/X cos 6); 


±j~rjh cos 6 






_—jh sin 6; Tj—jh cos 0; 


- 1 =F h (cos 6+ sin 2 8/X cos 


0). 






713 









(34) 



(35) 



(3«) 



with h defined as 

h=XY/2m 2 . (37) 

In the equations above, the upper sign is always associated with index 1 and the lower sign 
with the index 2. Here, however, the terms associated with the former index represent the 
ordinary wave when 0<9O°, but represent the extraordinary wave when 0>9O°. 

The exponents in (32) above are 

—jkN lt2 (R)'R=—jkpm(l±h cos e) (38) 

and the spherical polar components of the radiation vectors are calculated from those of the 
source currents through 

V It ,(fi) = {[{ i(R')e' kN ^ R) ' R 'p'* sin e'dp'dB'd*' (39a) 



and 



W li2 (£)= f(f M(R')e m ^ m '*'(>'* sin V&p'dB'ty (39b) 



with the exponent given by 

jkN h2 (R) • R f =jkp'm[sin Q' sin 6 cos (</ — p)+cos 0' cos d±h cos 6'] 

=jkp'm cos \l/±jkmhz' , (40) 

where \p is the angle between R and R' . 

Equations (32) through (40) constitute the desired approximate solution to the problem of 
calculating the radiation field of given source currents fixed in a lossless magneto-ionic medium 
and radiating at a frequency above the plasma frequency and well above the gyrofrequency. 

4. Application of the Solution to Radiation From Arrays 

The first order QL solution (32) holds over most of the region surrounding the source 
currents when the assumed conditions are met. It is interesting to consider the significance 
of this solution as it applies to the radiation from an array of similar sources. For brevity, 
only electric type sources are discussed. 

The current density in an array of n sources located at R x , R 2 , . . . R n can be written as 

n 

J(R)=J2.f i J s (R-R ( ) (41) 

i = l 

where J S (R) represents the current density of one of the sources, and the complex constants 
1i are feeding coefficients that describe the relative amplitudes and phases between the sources. 
Substituting into (39), the radiation vectors for such an array are 

Y 1>2 (R)=^fie jkp i mC0B *i± jkmhz i' J ( j s (B f )e ikp ' mC0B ^ jkmhz, p f2 sin B'dp'djd'db' . (42) 

The existence of the two radiation vectors V x and V 2 above indicates that the patterns of 
the two characteristic waves emitted by the array might, in general, just as well be treated 
separately. However, when the array is so limited in size or is so oriented that the term 
(kmhz) produces negligible change in phase over the array and its sources, then the radiation 
vectors V x and V 2 are equal, and the electric intensity (32a) can be written as 

E(JS)=^ e-^cos (kpmh cos 6) (A.+A^^-j sin (kpmh cos 6) (A 1 -A 2 )l2]\ (43) 

with 

y=^/«g**» fl0 "*< • ([{j s (R')e jk >' mc0S +p' 2 sin B'dp'dB'db' . (44) 
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The eq (44) for the radiation vector V above is identical with the formula that would give the 
radiation vector of the array it it were immersed in an isotropic medium having an index of 
refraction equal to m. It is equal to the product of the complex space factor of the array 
[Schelkunoff, 1943, sec. 9.3] times the radiation vector of the source. From this it follows that 
a principle of pattern multiplication [Kraus, 1950] can be stated for the array considered above 
as follows: 

In a lossless magneto-ionic medium at frequencies above the plasma frequency and well 
above the gyro frequency, the total field pattern of an array composed of similar sources and 
over which a phase change of magnitude (kmh) radians per unit length in the z direction can 
be neglected is equal to the product of the pattern of an array of point sources calculated by 
neglecting the effect of the static magnetic field and the pattern of an individual source in the 
anisotropic medium. The electric intensity from which the latter pattern is calculated is given 
by (43) using the radiation vector of an individual source in place of V. This electric intensity 
differs from that existing in an isotropic medium both because of the factors cos (Jcpmh cos 6) 
and sin (kpmh cos 6) appearing in (43) and because of the presence of terms containing the 

A A 

factor h in the matrices Ai and A 2 in (43). If the latter terms are small enough to be neglected, 
then the only effect of the anisotropy of the medium on the radiation field of the source and 
consequently of the entire array is the Faraday rotation due to the cosine and sine factors. 

In dealing with arrays in which the phase shift due to the terms (kmhz') in (40) is not 
negligible, the two radiation vectors V x and V 2 are distinct. Consequently, the radiation 
patterns of the associated characteristic waves cannot be combined as was done above. How- 
ever, a principle of pattern multiplication can still be stated for each characteristic wave. 
Since the term (kmhz') is independent of the direction of 7?, the phase shift due to this term can 
be associated with that of the feeding coefficients of the array making the following statement 
possible: 

Each characteristic wave emitted by an array of sources immersed in a lossless magneto- 
ionic medium and radiating at a frequency above the plasma frequency and well above the 
gyrofrequency has a pattern equal to the product of the pattern of the characteristic wave 
emitted by an individual source times a space factor of the array. The latter can be calculated 
by assuming the array to be immersed in the plasma with no static magnetic field present, 
provided that the feeding coefficients of the array are changed by a compensating phase shift 
of ±hnh radians per unit length in the z direction. 

5. Conclusions 

The electromagnetic field radiated by either an electric or a magnetic type source embedded 
in a lossless magneto-ionic medium has been given directly in terms of the components of two 
radiation vectors. One of the vectors gives the field due to the ordinary wave while the other 
gives the field due to the extraordinary wave. The components of the radiation vectors are 
obtained by an integration over the source. 

The approximations used to derive the solution restrict the application of the latter to 
cases where the source frequency is above the plasma frequency and is high compared to the 
gyromagnetic frequency of the medium. Also the relative geometrical placement between the 
source and the field point must be such that the QL approximation is valid. The Appleton- 
Hartree equations for a lossless magneto-ionic medium are used. 

An examination of the solution shows that the principal effect of the anisotropy of the 
medium on the radiation field of a source of limited extent is a Faraday rotation of the field 
that would exist if the medium were made isotropic by setting its static magnetic field equal 
to zero. There is in addition, however, some first order distortion of the pattern of the source 
proportional to the strength of the static magnetic field. The space factor of the pattern of 
an array of similar sources is affected only by the incremental Faraday rotation between the 
elements of the array. When the extent of the array is such that this incremental rotation is 
negligible, the space factor of the array is identical with that obtained by treating the medium 
as isotropic. When the rotation is not negligible a distinct array space factor is obtained for 
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each of the characteristic waves. These space factors can also be calculated, however, from 
the formulas for the pattern of the array in an isotropic medium by making compensating 
changes in the phases of the feeding coefficients of the array elements. 

6. References 

L. G. Abraham, Jr. (Aug. 1, 1953), Extensions of the magneto-ionic theory for radio-wave propagation in the 
ionosphere including antenna radiation and plane wave scattering, Cornell Univ. School of E. E. Tech. 
Kept. No. 13 on studies on propagation in the ionosphere. 

M. P. Bachynski (Dec. 1961), Plasma physics — an elementary review, Proc. IRE 4, 1751-1766. 

F. V. Bunkin (Feb. 1957), On radiation in anisotropic media, Soviet Phys .— JETP 32, 338-346. 

Y. Chow (July 1962), A note on radiation in a gyro-electric-magnetic medium — an extension of Bunkin's cal- 
culation, IRE Trans. Ant. Prop. AP-10, 464-469. 

G. W. Ford (1961), Electromagnetic radiation from a source in a plasma, Ann. Phys. 16, 185-200. 

H. Kogelnik (Sept.-Oct. 1959), On electromagnetic radiation in magneto-ionic media, J. Res. NBS 64D (Radio 

Prop.), 515-522. 
J. D. Kraus (1950), Antennas, sec. 4-3 (McGraw-Hill Book Co., Inc., New York, N.Y.). 
H. H. Kuehl (Sept. 1962), Electromagnetic radiation from an electric dipole in a cold anisotropic plasma, Phys. 

Fluids 5, 1095-1103. 
M. J. Lighthill (1960), Studies on magneto-hydrodynamic waves and other anisotropic wave motions, Phil. 

Trans. Roy. Soc. A252, 397-430. 
J. W. Marini (1963), Radiation from sources in magneto-ionic media, Doctoral Thesis, University of Maryland. 

Obtainable from University Microfilms, Inc., 313 N. First St., Ann Arbor, Mich. 
J. A. Ratcliffe (1959), The magneto-ionic theory and its applications to the ionosphere (Cambridge Univ. Press, 

New York, N.Y., and London, England). 
S. A. Schelkunoff (Oct. 1939), A general radiation formula, Proc. IRE 27, 660-666. 

S. A. Schelkunoff (1943), Electromagnetic waves, sections 4.1 and 4.4 (D. Van Nostrand Co., New York, N.Y.). 
R. Shore and G. Meltz (Jan. 1962), Anisotropic plasma-covered magnetic line source, IRE Trans.- Ant. Prop. 

AP10, 78-82. 
H. Suhl and L. R. Walker (May 1954), Topics in guided wave propagation through gyromagnetic media, Bell 

Sys. Tech. J. 33, 579-659. 

Additional Reference 

S. Ramo and J. R. Whinnery (1959), Fields and Waves in Modern Radio, Ch. 12 (John Wiley & Sons, Inc., 
New York, N.Y.). 

(Paper 67D6-300) 



716 



